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Abstract. In previous work, the first author developed an algorithm for the computation 
of Hilbert modular forms. In this paper, we extend this to all totally real number fields of 
even degree and nontrivial class group. Using the algorithm over Q(VTO) and Q(v°5) and 
their Hilbert class fields, we present some new instances of the conjectural Eichler-Shimura 
construction for totally real fields, and in particular find new examples of modular abelian 
varieties with everywhere good reduction. 



Introduction 

Let F be a totally real number field of even degree. Let B be the unique quaternion al- 
gebra over F which is ramified exactly at all infinite places. By the Jacquet-Langlands cor- 
respondence (Jacquet and Langlands [H Chap. XVI] and Gelbart |6J), computing Hilbert 
modular forms over F amounts to computing automorphic forms on B. In [2] and [3], 
the first author presented an algorithm that exploits this correspondence. The algorithm 
adopts an alternative approach to the theory of Brandt matrices on B that is compu- 
tationally more efficient than the classical one. Both papers considered only fields with 
narrow class number one. One technical difficulty arising from nontrivial class groups is 
that ideals in B are no longer free Op-modules. This is now handled smoothly in the 
package for quaternion algebras over number fields included in the Magma computational 
algebra system pQ (version 2.14). Our computations rely heavily on this package, in which 
algorithms from [15] and [9] are implemented. There are not many explicit examples of 
Hilbert modular forms in the nontrivial class group case in the literature. Okada [11] pro- 
vides several examples of such forms of level 1 and parallel weight 2 on the quadratic fields 
Q(v / 257) and Q(V40l), computed using explicit trace formulae. However, it would be 
difficult to extend this approach to arbitrary totally real number fields, levels and weights. 
In this paper we present a general algorithm that is practical for a large range of fields 
and levels. This opens the possibility of experimenting systematically, especially over fields 
with nontrivial class group, and we hope this will shed new light on the theory of these 
objects. 

The paper is organized as follows. Section [1] contains the necessary theoretical back- 
ground. In section [2] we state the general algorithm, and describe some improvements to 
its implementation. Section [3] provides some numerical data over the real quadratic fields 
Q(\/l0) and Q(y / 85) and their Hilbert class fields. In section H] we use this data to give 
new examples of the Eichler-Shimura construction over totally real number fields. 
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1. Theoretical background 

Our aim is to compute spaces of Hilbert modular forms (as Hecke modules). By the 
Jacquet-Langlands correspondence, this is equivalent to computing spaces of automorphic 
forms Mj?(N) on some suitable quaternion algebra B (as Hecke modules). In this section 
we define Mf(N), and then explain how to compute it in terms of spaces M^(R, N) of 
automorphic forms on quaternion orders R C B. A good reference for the material on 
Hilbert modular forms is Taylor [TJ]. For the theory of Brandt matrices, we refer to 
Eichler [5], and also to [3] (which translates it into the adelic framework used here). 

We fix a totally real number field F of even degree. We let I be the set of all real 
embeddings of F; and for any r G J, we denote the corresponding real embedding by 
a I— > a T , a G F. Also, we let Op be the ring of integers of F, and fix an integral ideal 
iV of F. We let B be the unique (up to isomorphism) quaternion algebra ramified at 
the infinite places of F. We fix a maximal order R of B. We fix a Galois extension K 
of F contained in C, which splits B. We also fix an isomorphism B <S>f K = M^A) 7 , 
and let j : B x GL^C) 7 be the resulting embedding. For each prime p of F, we 
fix a local isomorphism B p = M^-Fp) such that R p = M^CV.p)- Combining these local 
isomorphisms, one obtains an isomorphism B = M 2 (-F) under which R = M 2 (C^), where 
F and Of are the finite adeles of F and Of respectively. We define the compact open 
subgroup Uq(N) of R x by 

^ := {(c d) eGL 2 (a F ): c = 0(JV)}. 

Let CIb denote a complete set of representatives of all the right ideal classes of R. Then 
CIb is in bijection with B X \B X /R X . We choose a finite set of primes S that generate the 
narrow class group Clp and such that q \ N for any q G S. For any a G CIb, we let R a be 
the left (maximal) order of Q. Applying the strong approximation theorem, we may choose 
the representatives d e CIb so that the prime divisors of Na lie in S. Then there are 
well-defined surjective reduction maps R x — ► GL 2 (Cf/^), that all differ by conjugation 
in G~L2(Of/N). This gives a transitive action of each R x on P 1 (Of/N). 

Fix a vector k G Z 7 such that k T > 2 for all r, with all the components having the same 
parity. Set t = (1, . . . , 1) and m = k — 2t, then choose v&Z 1 such that each v T > 0, v T = 
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for some r, and m + 2v = fit for some non-negative /16Z. For every non-negative integer 
a, b G Z, we let S a; fc(C) denote the right M 2 (C)-module Sym a (C 2 ) (the a th symmetric 
power of the standard right M2(C)-module C 2 ) with the M2(C)-action: 

x ■ m :— (det m) b xSym a (m) . 

Then, we define the weight representation by 

^fc = ®S mr , Vr (C). 

rel 

The space of automorphic forms of level N and weight k on B is defined as 
Mf (AT) := {/ : B x /U (N) - L k _: = f \ 7 G B x } , 

where /||fc7(x) := f(jx)~/. 

Remark 1. By the Jacquet-Langlands correspondence [7J Chap. XVI], there is an isomor- 
phism of Hecke modules between M®(N) and M k (N), the space of Hilbert modular forms 
of weight k and level N over F. Therefore our task is now to compute M®(N) as a Hecke 
module. 

We define the space of automorphic forms of level N and weight k on the order R n by 
Mfc(R a , N):={f: P\O f /N) - L k : f\\ kl = /, 7 e K } ■ 

For each a, b G C7b and any prime p G F, put 

B^p, a, b) := R a *\ G ab- 1 : ( ^ (m)) i = p\ , 

[ nr(a)nr(b) 1 J 

where R x acts by multiplication to the left. We define the linear map 

T a , b (p): M K (R b ,N) -> Mfc(R fl ,JV) 

«ee( s )(p, a, 6) 

Proposition 1. There is an isomorphism of Hecke modules 

K{N) -> Mfc(R a , iV), 

aeCZ B 

where the action of the Hecke operator T(p) on the right hand side is given by the collection 
of linear maps (T a> &(p)) ; 0, b G CI-q. 

Proof. The Brandt matrices for a totally real number field are defined in Eichler j5j, using 
the global language. In that language, it is not hard to see that the action of the Hecke 
operator T(p) on the space M®(1) is determined by the collection of sets 6^ 5 ^(p, a, b), 
a, b G CIb, for each prime p. Let G = ReSi?/Q(B x ) and apply [U Theorem 1] to obtain 
Proposition [TJ 
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Alternatively, one can observe that the proof of \3, Theorem 2] does not use the restriction 
on the class number of F. And so, Proposition [T] is simply a translation of 0, Theorem 2] 
in global terms. □ 



2. Algorithmic issues 

In the case of real quadratic fields, our algorithm has been discussed in J2j sec. 2] and [H 
sec. 6]. In this section we state the algorithm in full generality, for any totally real number 
field F of even degree and any weight and level. We then discuss some implementation 
issues. 

As above, B will be the quaternion algebra ramified at the infinite places of F. Our 
goal is to compute the space Mf*(N) of automorphic forms on B of weight k and level N, 
where N is an integral ideal in F. First choose a (reasonable) bound b G N, and a set of 
prime ideals S not dividing N that generate Clp. The more expensive tasks are done in 
an initial precomputation, since these depend only on F and not on the weight or level. 

Precomputation. 

(1) Compute a maximal order R of B. 

(2) Compute a complete set CIb of representatives a for the right ideal classes of R, 
chosen so that all prime factors of Na lie in S. 

(3) For each representative a G CIb, compute its left order R a , and compute the unit 
group r a of R a . 

(4) Compute the sets 9^(p, a, b), for all primes p with Np < b and all a, b G CIb- 
(See Section I2TT1 for details.) 

Algorithm. 

(1) Compute splitting isomorphisms R* = GL 2 (0f iP ), for each prime p | N. 

(2) For each a G CIb, compute M k (R a , N) as a module of coinvariants 

il4(R Q , N) = K[V\O f /N)} ® L k /(x - 7 i,7e T a ). 

(3) Compute the direct sum 

Mf{N) = A4(R a , N). 

a&Cl B 

(4) For each prime p of F with Np < b, compute the families of linear maps (T 0i (,(p)). 
(These determine the Hecke operator T(p) as a block matrix.) 

(5) Find a common basis of eigenvectors of M®(N) for the T(p). 

(6) If Step (5) does not completely diagonalize M^(N), increase b and extend the 
precomputation, obtaining 6^(p, a, b) for Np < b. Then return to Step (4). 

Remark 2. In practice, it is extremely rare that one resorts to Step (6) since very few 
Hecke operators T(p) are required to diagonalize the space M®(N). In the cases we tested, 
which included levels with norm as large as 5000, we never needed more than 10 primes. 
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The steps in the main algorithm involve only local computations and linear algebra, 
whereas several steps in the precomputation involve lattice enumeration. For a given field 
F, the precomputed data G^^p, a, b) can be re-used for all levels that are coprime to 
the primes in S, and all weights. If one wishes to compute for all levels up to some large 
bound, one may simply choose the primes in S to be larger than the bound. 

2.1. Computing 9^(p, a, b). Recall 9^(p, a, b) C R^ab" 1 (as defined in Section [TJ). 

Lemma 2. The correspondence u <-» u~ x a gives a bijection between 0^(p, a, b) and the 
set of fractional right K-ideals c > b such that nr(b) = nr(c)p and c = a as right K-ideals. 

Proof. Consider all fractional right R-ideals c that are isomorphic to o; these are precisely 
the ideals u~ l a for u G B x . Note that u~ x a = v~ l a if and only if v G R*m. It is clear 
that u~ l a contains b if and only if u G ab -1 , and that nr(b) = nr(u _1 o)p if and only if 
(nr(it)) = nr(a)nr(b) _1 p. The lemma follows. □ 

Algorithm. This computes B^(p, a, b) for all o G Cl&, where p and b are fixed. 

(1) Compute the fractional right R-ideals c > b with nr(b) = nr(c)p. (The number of 
these is Np + 1.) 

(2) For each such c, compute the representative a G CI-q and some u G B such that 
c = n _1 a. Append u to 9^(p, a, b). 

Remark 3. One sees that for each fixed p and b, 

#© (5) (P, a, b) = Np + l, 

aeC7 B 

however this fact is not used in the algorithm. 

Step (1) is a local computation; the ideals are obtained by pulling back local ideals under 
a splitting homomorphism R p = M^-Fp). Step (2) is the standard problem of isomorphism 
testing for right ideals, and we discuss an improvement to the standard algorithm for this 
in the next section; the complexity of each isomorphism test will not depend on p. 

2.2. Lattice-based algorithms for definite quaternion algebras. In this section, we 
let B be any definite quaternion algebra over a totally real number field F, and let R be a 
(maximal or Eichler) order of B. Two basic algorithmic problems are: 

(1) to find an isomorphism between given right R-ideals a and b, and 

(2) to compute the unit group of R (modulo the unit group of Of)- 

The standard approach to both problems (as in [15]) reduces them to the following. 

General Problem: Let L = Z" be a lattice contained in B (not necessarily of full rank). 
Given a totally positive element a G F , compute all x G L with nr(x) = a. 

For isomorphism testing, L is the fractional ideal ab -1 and a is a generator of nr(ob _1 ). 
(It suffices to consider a finite set of possibilities for a.) For computing units, L is R (or 
possibly an CV-sub module of R known to contain a unit), and a is a unit of Of- 
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One may solve the general problem by considering the positive definite quadratic form 
on L given by Tr(nr(rr)). (Note that its values are positive since nr(x) is a totally positive 
element of F, for all 7^ x G B). We capture all x G L with nr(x) = a by enumerating 
all x for which the quadratic form takes value Tr(a), using standard lattice algorithms. 
The drawback is that Tr(ai) might not be particularly small in relation to the determinant 
of the lattice (even when a is a unit), in which case the lattice enumeration can be very 
time-consuming. 

We now present a variation which avoids this bottleneck. For any nonzero c G F, we 
may instead consider the lattice cL C B, again under the positive definite quadratic form 
given by Tr(nr(x)). We clearly capture all x G L with nr(a;) = a by enumerating all y G cL 
with Tr(nr(y)) = Tr(c 2 a) and taking x = y/c. In the special case that c G Q, this merely 
rescales the enumeration problem. However, we will see that c G F may be chosen so that, 
in applications (1) and (2) above, one is looking for relatively short vectors in the lattice. 

Note that det(cL) = |N(c)| dim ( L )/ dcgF det(L). Heuristically, as c varies, the complexity 
of the enumeration process will be roughly proportional to the number of lattice elements 
with length up to the desired length. Asymptotically this number equals 

Tr(c 2 a )dim(L)/2 Tr(c2a) dim(L)/2 Tf ^dimQL)^ N ( a )dim(L)/2 dcg F 

det (cL) : : \N{c)\ dim ( L V dc z F det{L) ' : N(c 2 a) dim ( L V 2dc § F det(L) ' 

Given that a is totally positive, Tr(c 2 a)/N(c 2 o;) 1//degF cannot be less than degF, and is 
close to degF when all the real embeddings of c 2 a lie close together. It is straightforward 
to find c G Of with this property, as follows. First fix a Z-module basis bas(Op) of Of- 

Algorithm. Choose a large constant C . 

(1) Calculate r; L := C '/ ^^/a^Ja) (note that the real embeddings of a are positive). 

(2) Represent the vector (r*) in terms of the basis bas(Op), then round the coordinates 
to integers, thus obtaining an element c G Of- 

Lemma 3. Given any totally positive element a G F, and any e > 0, we can find c G F 
such that Tr(c 2 a)/N(c 2 a) 1 / dc §^ < degF + e. 

Proof. In the notation of the algorithm, we fix a and let C — > 00, regarding r\ G K and 
c G Op as functions of C. Since we use a fixed basis of Op, &i{c) — r\ is bounded by a 
constant independent of C. Therefore as C — > 00, a,i(c 2 a) = rfa^a) + 0(C) = C 2 + 0(C) . 
This implies that for any i and j, the ratio ai(c 2 a)/aj(c 2 a) — > 1 as C — > 00, and the lemma 
follows. □ 

The complexity of the enumeration thus depends on the ratio N(a) dim( ^ L - )/ ' 2degF /det(L). 
In both the applications above, this ratio is small: in computing units, a is a unit, and in 
isomorphism testing, a generates the fractional ideal nr(L) where L = ab^ 1 . 

3. Examples of Hilbert modular forms 

In this section we give some examples computed using our algorithm, which we have 
implemented in Magma (and which will be available in a future version of Magma). 
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N(p) 


P 


ft 


h 




/s 


u 


3 


(3, 2uj 85 ) 


4 


-4 


2\ 


/-[ 




3 


(3,4 + 2^ 85 ) 


4 


-4 


-2\ 


/-L 


P 


4 


(2,0) 


5 


5 




1 


-/? 3 + 3 


5 


(5, -1 + 2uj 85 ) 


6 


-6 







-/? 3 + 4/3 


7 


(7, 2lu 85 ) 


8 


-8 


-2\ 


f—Y 


/3 3 — 5/3 


7 


(7, 12 + 2cu 85 ) 


8 


-8 


2v 


f—Y 


/3 3 -5/3 


17 


(17,-l + 2^ 85 ) 


18 


-18 







2/3 3 - 14/3 


19 


(19,2 + 2w 85 ) 


20 


20 




-4 


2 


19 


(19, 15 + 2oo 85 ) 


20 


20 




-4 


2 



Table 1. Hilbert modular forms of level 1 and parallel weight 2 over 
Q(v / 85)- (The minimal polynomial of (3 is x 4 — 6x 2 + 2). 



3.1. The quadratic field Q(v / 85). Let F be the real quadratic field Q(a/85). The class 
number of F is the same as its narrow class number: hp = hp = 2. The maximal order in 
F is Op = Z[cu 85 ], where lu 85 = Let B/F be the Hamilton quaternion algebra, i.e., 

the F-algebra given by 

B := F © Fi © Fj © Fk, with i 2 = -1, j 2 = -1 and fc = ij. 

Since the prime 2 is inert in F, the algebra B is ramified only at the two infinite places. 
Using Magma, we find that the class number of B is 8. The Hecke module of Hilbert 
modular forms of level 1 and weight (2, 2) over F is therefore an 8-dimensional Q-space, 
and it can be diagonalized by using the Hecke operator T2. There are two Eisenstein series 
and two Galois conjugacy classes of newforms. The eigenvalues of the Hecke operators for 
the first few primes are given in Table [T] (only one eigenform in each Galois conjugacy class 
of newforms is listed). 

The Hilbert class field of Q( v / 85) is H := Q{V5, VTf) = Q(a), where a 4 - 4a 3 - 5a 2 + 
18a — 1 = 0. The narrow class number of H is 1, and B ®p H (the quaternion algebra 
over H ramified at the four infinite places) has class number 4. Thus the space of Hilbert 
modular forms of level 1 and weight (2, 2) is 4-dimensional. The eigenvalues of the Hecke 
action for the first few primes are listed in Table El There is one Eisenstein series and two 
classes of newforms. Elements of Oh are expressed in terms of the integral basis 

1, -(a 3 -3a 2 -5a + 10), -(-a 3 + 3a 2 + 11a - 10), -(-a 3 + 14a + 5), 
6 6 6 

which we use to write generators of the ideals in the table. 

We also computed some spaces over Q(y / 85) with nontrivial level. The dimensions of 
the spaces with prime level of norm less than 100 are given in Table [3j (It suffices to 
consider just one prime in each pair of conjugate primes, and for the precomputation we 
took S = (3, —1 + cl> 85 ).) For example, for level p = (5, a/85), M 2 (p) has dimension 20, and 
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N(p) 


P 


h 


h 


h 


4 


[1,-1,0,1] 


5 


1 


3 + F 


4 


[0,2,-1,1] 


5 


1 


s + P' 


9 


[0,1,-1,0] 


10 


2 


F 


9 


[1,-1,-1,0] 


10 


2 


F 


19 


[0,1,0,-1] 


20 


-4 


2 


19 


[-1,2,0,1] 


20 


-4 


2 


19 


[1,-1,-1,1] 


20 


-4 


2 


19 


[-1,2,-1,1] 


20 


-4 


2 



Table 2. Hilbert modular forms of level 1 and parallel weight 2 over the 
Hilbert class field H of Q(y / 85)- (The minimal polynomial of f3' is x 2 +6x+2). 



N(p) 


dimM 2 (p) 


dimS 2 (p) 


dim S% ew (p) 


3 


16 


14 


8 


4 


24 


22 


16 


5 


20 


18 


12 


7 


32 


30 


24 


17 


56 


54 


48 


19 


68 


66 


60 


23 


72 


70 


64 


37 


124 


122 


116 


59 


180 


178 


172 


73 


232 


230 


224 


89 


272 


270 


264 


97 


304 


302 


296 



Table 3. Dimensions of spaces of Hilbert modular forms over Q(a/85) with 
weight (2, 2) and prime level of norm less than 100 



the Hecke operator T q with q = (7, 2cj 85 ) acting on M 2 (p) has characteristic polynomial 

(x -8)(x + 8)(x 2 + 4) 2 (a; 4 - 10x 2 + 18) 2 (x 6 + 28x 4 + 104a; 2 + 100). 

Comparing this with the space M 2 ((l)) of level 1, on which T q has characteristic polynomial 

(x - 8)(x + 8)(x 2 + 4)(x 4 - 10x 2 + 18), 

one sees that the Hecke action on the subspace of newforms M 2 (p) is irreducible, and the 
cuspidal oldform space embeds in M 2 (p) under two degeneracy maps (as expected). 

3.2. The quadratic field Q(v / 10). Let F be the real quadratic field Q(v / 10)- The Hilbert 
class field of F is H := Q(v / 2, V§) — Q(«), where the minimal polynomial of a is x A — 
2a; 3 — 5a; 2 + 6a; — 1. The narrow class number of H is 1. We computed the space of Hilbert 
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N(p) 


P 


p 

h 


p 

h 


p 

h 


2 


(2,cu 40 ) 


-3 


3 




3 


(3, ^40 + 4) 


-4 


4 


V2 


3 


(3, ^40 + 2) 


-4 


4 


V2 


5 


(5, u) 4 o) 


-6 


6 


-2V2 


13 


(13, lj 40 + 6) 


-14 


14 





13 


(13, uj 40 + 7) 


-14 


14 





31 


(31,w 40 + 14) 


32 


32 


4 


31 


(31,^40 + 17) 


32 


32 


4 


37 


(37,^40 + H) 


-38 


38 


6V2 


37 


(37, oo m + 26) 


-38 


38 


6V2 



Table 4. Hilbert modular forms of level 1 and weight (2, 2) over Q(V4~0). 



N(P) 


P 


fi 


h 


4 


[0,0,1,0] 


5 


-2 


9 


[1,1,-1,0] 


10 


-4 


9 


[0,1,-1,1] 


10 


-4 


25 


[1,-2,0,0] 


26 


-2 


31 


[1,1,1,-1] 


32 


4 


31 


[1,-1,-1,-1] 


32 


4 


31 


[1,1,-1,1] 


32 


4 


31 


[-3,2,-1,0] 


32 


4 



Table 5. Hilbert modular forms of level 1 and weight (2, 2) over the Hilbert 
class field H of Q(v / 10). 

modular forms of level 1 and weight (2, 2) over F and H, and the Hecke eigenvalues for 
the first few primes are listed in Table H] and Table [5] (only one eigenform in each Galois 
conjugacy class of newforms is listed). Elements of Oh are expressed in terms of the 
integral basis 

1, -{2a 3 - 3a 2 - 10a + 7), -{-2a 3 + 3a 2 + 13a - 7), -{-a 3 + 3a 2 + 5a - 8). 
3 3 3 

4. Examples of the Eichler-Shimura construction 

In the study of Hilbert modular forms, the following conjecture is well-known. We refer 
to Shimura [12] or Knapp [10] for the classical case, and to Zhang [16] and references 
therein for the number field case. 

Conjecture 4 (Eichler-Shimura). Let f be a Hilbert newform of level N and parallel weight 
2 over a totally real field F. Let Kf be the number field generated by the Fourier coefficients 
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of f . Then there exists an abelian variety Af defined over F such that Kf 
and 

L(A f , s)= J] L(f , s), 

cr£Gai(K f /Q) 

where f a is obtained by letting a act on the Fourier coefficients of f . 

In the classical setting, namely when F = Q, this is a theorem known as the Eichler- 
Shimura construction. In general, many cases of the conjecture are also known. (See, for 
example, Zhang [16J and references therein). In those cases the abelian variety Af is often 
constructed as a quotient of the Jacobian of some Shimura curve of level N. The case 
where [F : Q] is even and Af has everywhere good reduction poses a different challenge. 
In this section, we provide new examples of such Af. 

Remark 4. We refer back to the final paragraph of section 13.11 The characteristic poly- 
nomials given there, viewed in terms of Conjecture HI indicate that the newform part of 
M 2 (p) corresponds to a simple abelian variety of dimension 6. 

4.1. The quadratic field Q(v / 85). Keeping the notation of subsection 13.11 let E/H be 
the elliptic curve with the following coefficients: 





a 2 a 3 


04 


a 6 


E: [1,0,0,1] 


[0,-1,0,-1] [0,1,1,0] 


-5,-6,-1,0] 


-8,-7,-3,2] 



This curve has everywhere good reduction. The restriction of scalars A = lies h/f{E) is an 
abelian surface over F, also with everywhere good reduction. 

Remark 5. The j-invariant of E is 64047678245 - 125343498 15lu 85 e F, and in fact E is 
./^-isomorphic to its conjugate under Gal(H/F). Therefore A is if-isomorphic to E © E. 
Let E' denote one of the other two conjugates, which have j-invariant 51513328430 + 
12534349815u;85; there is an isogeny of degree 2 from E to E' . The restriction of scalars 
ResH/F(E') over F also has good reduction everywhere; it is if-isomorphic to E' ® E' , and 
is therefore isogenous to A. 

To establish the modularity of E and A, we will apply the following result of Skinner 
and Wiles. Here we state the nearly ordinary assumption (Condition (iv)) in a slightly 
different way. 

Theorem 5. [T3| Theorem A] Let F be a totally real abelian extension o/Q. Suppose that 
p > 3 is prime, and let p : Gal(F/F) — > GL 2 (Q P ) be a continuous, absolutely irreducible 
and totally odd representation unramified away from a finite set of places of F. Suppose 
that the reduction of p is of the form p ss — xi © X2, where Xi ^nd X2 o^re characters, and 
suppose that: 

(i) the splitting field F{x\ IX2) of X1/X2 is abelian over Q, 
(H) (xi/x*)\d v 1 f or e «c/i v I p, 



End(A/)®Q 
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(Hi) p\i v = I p ^ J /or eac/i prime v \ p, 

(iv) det p = V e p "S wi/i k > 2 an integer, ip a character of finite order, and e p the 
p-adic cyclotomic character. 

Then p comes from a Hilbert modular from. 

Proposition 6. a) The elliptic curves E is modular and corresponds to / 2 in Tabled 

b) The abelian surface A is modular and corresponds to f% in TableUi 

Proof, a) Let pe,3 be the 3-adic representation attached to E, and pE,3 the corresponding 
residual representation. Also, let p be any prime above 3. Using Magma, we compute the 
torsion subgroup E(H) tors = Z/2 ©Z/2, and the trace of Frobenius a p (E) = 2. The latter 
implies that the representation pe,3 is ordinary at p. By direct calculation, we find that 
j{E) is the image of a if -rational point on the modular curve X (3): 

(t + 27) (r + 3) 3 
j(E) = ± ^ '-, where r = [2166, 527, -527, 1054]. 

T 

This implies that E has a Galois-stable subgroup of order 3, so the representation Pe,3 
is reducible. Since it is ordinary, there exist characters \i X 1 unramified away from p | 3, 
with x unramified at p, such that p s £ 3 = \ ® x' an d XX! — e 3 is the mod 3 cyclotomic 
character. The field H(x/x') is clearly abelian. Therefore the representation 3 satisfies 
the conditions of Skinner and Wiles, and E is modular. Comparing traces of Frobenius 
with the eigenvalues given in Tabled we see that the corresponding form is /g. 

6) Let BC(/3) be the base change from F to H of the newform f 3 in Tabled) Since the 
Hilbert class field extension H/F is totally unramified, the form BC(/3) has level 1 and 
trivial character. By comparing the Fourier coefficients at the split primes above 19, we 
see that BC(/s) = f'2 in Tabled The result then follows from properties of restriction of 
scalars and base change. □ 

Remark 6. To find E, we reasoned as follows. The eigenvalues of fi in Table [2] suggest 
that the corresponding curve admits a 2-isogeny. The curve must have good reduction 
everywhere, and so must its conjugates; if these are also modular, then they share the 
same L-series and are therefore isogenous to each other. This would mean the curve comes 
from an if -rational point on Xo(2) whose j-invariant is integral. Using a parametrisation 
of X (2), we searched for such points. We would like to thank Noam Elkies for suggesting 
this approach. 

Remark 7. If we assume Conjecture HI then there exists a modular abelian surface A over 
H with real multiplication by Q(Vi) which corresponds to the form / 3 in Table [21 The 
restriction of scalars of A from H to F is a modular abelian fourfold with real multiplication 
by Q(/3) which corresponds to the form / 4 in Table [H 



12 LASSINA DEMBELE, STEVE DONNELLY 

4.2. The quadratic field Q(V^IO). Keeping the notation of subsection 13.21 let E/H be 
the elliptic curve with the following coefficients: 



CLi 


a 2 a 3 


04 


a 6 


E : [0,0,1,0 


1,0,1,-1] [0,1,0,0 


-15,-44,-21,-26 


-91,-123,-48,-97] 



This is an elliptic curve with everywhere good reduction over H. In contrast with the 
previous example, the four Galois conjugates have distinct j-invariants. The restriction of 
scalars A = Res#/i?(i?) is an abelian surface over F with everywhere good reduction. 

Proposition 7. The elliptic curve E/H and the abelian surface A/F are modular; E 
corresponds to fa in Tabled and A corresponds to fa in Table^ 

Proof. Let pe,3 be the 3-adic representation attached to E, and pe,3 its reduction modulo 
3. Then, p E ,3 is reducible since 

= (r + 27 ^ T + 3)3 , where r = [5, 52, -18, -26]. 

As before, it is easy to see that pe,z satisfies the conditions of Skinner and Wiles. So E 
is modular, and hence A is also modular. Comparing traces of Frobenius with Fourier 
coefficients, it is easy to see which forms in the tables they correspond to. 

Alternatively, we could consider the 7-adic representation Pej- Its reduction mod 7 is 
reducible since the point ([16,23,9,18] : [-157,-268,-119,-184] : [1,0,0,0]) is an ir- 
rational point of order 7 on E. Furthermore, for any prime p | 7, we have a p (E) = 8, and 
it is easy to see that p E j satisfies the conditions of Skinner and Wiles. □ 

Remark 8. It was shown by Kagawa [HI Theorem 3.2] that there is no elliptic curve with 
everywhere good reduction over Q(vTo). Our results show that if we assume modularity 
in addition, there is only one such simple abelian variety: an abelian surface with real 
multiplication by Z[\/2]- 

Remark 9. To find E, we were again assisted by the eigenvalues of the corresponding form 
fa in Table O which suggest that E has an ^/-rational point of order 14. The modular 
curve A (14)/Q is an elliptic curve (14A1 in Cremona's table), which (using Magma) was 
found to have rank 1 over H and also rank 1 over Q(\/l0); this enabled us to obtain a 
point of infinite order simply by finding a Q-rational point on the quadratic twist by y/l0. 
We considered curves corresponding to points of small height in X (14)(H), and twists of 
these curves, until we found one with good reduction everywhere. 
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